Scaling in the correlation energies of two-dimensional artificial atoms 
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We report an unexpected scaling in the correlation energy, E corr , of artificial atoms. Extensive 
numerical calculations of the ground-state energy of charged 2D-quantum dots show that the cor- 
relation energies scale as E corr /{Tiuj) « N 3 ^ A f corr {N 1 ^ i fi), where N is the number of electrons, the 
coupling constant /3 is the ratio between Coulomb and oscillator (hu)) characteristic energies, and 
fcorr is a universal function. An analytic expression for f CO rr is provided based on a two-parameter 
fit. In addition, analytic expressions for the correlation energy per particle and for the fraction of 
the total energy associated to the correlation energy are also provided. 



I. INTRODUCTION 

Artificial atoms - quantum dots (QD) [l], 0], i.e., 
nanoscopic semiconductor structures where a set of elec- 
trons is confined, offer wider possibilities of engineering 
their properties than real atoms. The atom size, for ex- 
ample, can be changed from a few nanometers to hun- 
dreds of nanometers by properly choosing the parame- 
ters. The degree of correlation of the electronic motion 
is also a property that can be widely modified because, 
in addition to Coulomb interaction, we can play with the 
confining potential. 

In the present paper, we study a very simplified quasi- 
twodimensional, isotropic, harmonic oscillator model of 
artificial atom. In spite of its simplicity, it is used very 
often to model real QD's The model can be charac- 
terized by two parameters, the number of electrons, N, 
and the ratio of Coulomb to harmonic typical energies 
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where e is the electron charge, m its effective mass in the 
semiconductor material, and e the dielectric constant. 

In a previous work Q , we showed that the total energy 
of such an artificial atom with dozens of electrons obeys 
the scaling relation dictated by Thomas-Fermi (TF) the- 
ory 



N 3/2 f gs (z). 
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The variable z = N 1 ^ 4 ^, combines in a particular way 
the number of electrons, N, and the coupling constant /3. 
The function f g8 is said to be universal in the sense that 
it depends on z. 

This result, in some sense, is not surprising, because 
TF is a mean-field theory, well describing properties of 
relatively "large" electron systems. It is amazing, how- 
ever, that the scaling holds in a wide range of j3, from 



the strong confinement limit (weak correlations, j3 — > 0) 
to the weak confinement regime (strong correlations, 
/? — > co, the so-called Wigner phase). 

In the present paper, we move a step further and check 
whether a scaling a la Thomas-Fermi, with different ex- 
ponents, holds also for the correlation energy 0-0], E corr , 
i.e., the difference between the total ground-state, E gs , 
and the Hartrec-Fock, Ehf, energies: 
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The fact that E corr scales is completely unexpected be- 
cause, by definition, correlation is precisely beyond mean 
field properties 0, B, 0-fH ■ 

We perform extensive calculations of the ground state 
of QD's with 6 < N < 56, and we find that 
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where f corr is a universal function. In addition, we find 
that the fraction of the total energy associated to E corr 
also scales in a universal way. 

The paper is organized as follows. In the next sec- 
tion we briefly summarize the computational methods 
employed. In Section IlIIl we present and discuss results 
from numerical calculations of E corr . Finally, concluding 
remarks are given in Section [IVI 



II. COMPUTATIONAL METHODS 

Extensive numerical calculations for charged quantum 
dots were performed. We follow standard procedures, 
which include variatonal Monte Carlo (VMC), density 
funtional theory (DFT), in particular the local den- 
sity approximation version (LDA), and both full (FCI) 
and truncated (2p2h CI) configuration interaction ap- 
proaches. 
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The LDA scheme is employed mainly to verify the per- 
formance of VMC in systems with large particle numbers, 
whereas the FCI (exact diagonalization) scheme is em- 
ployed, with the same objective, in the cases of small 
particle numbers. Additionally, we use a truncated 2p2h 
CI scheme in order to determine how much correlation is 
captured by it. The generalities of the implementations 
are briefly explained below. 



A. Hartree-Fock Scheme 

In our calculations we will represent the HF orbitals in 
the basis of Fock-Darwin states [l[ , 

fair) = ^2C al (j>i(f)xi, (5) 

i 

where Xi are the spin functions. The expansion coeffi- 
cients C a i and the energy e a can be obtained by solving 
the following eigenvalue problem (HF equations), 

^ y hjj C a j £aC a i 7 (6) 
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where the are the matrix elements of the HF self- 
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consistent Hamiltonian in the chosen basis, 

h* F = e i S ij + Ant E E[^ w l^> 

- {hul-lvJ^C^C^. (7) 

In Eq. §7$ the first term is diagonal, q denotes the 
single-particle oscillator energy e; = huj(2ki + \l\ + 1), 
and the second term accounts for direct and exchange 
Coulomb interaction between the electrons in the QD. 
The index p runs over the electron occupied HF orbitals 
and £/ denotes the Fermi level in the conduction band 
and p in t = (hu})/3. 

The HF equations, Eq. ©, are solved iteratively. 
Twenty oscillator shells (420 oscillator states) are used 
in the calculations. 



B. 2p2h Configuration Interaction Scheme 

In the truncated 2p2h configuration interaction 
scheme, the starting point Q is the Hartree-Fock solu- 
tion of the problem. Then, a basis of functions made up 
from (i) the Hartree-Fock state, \HF), (ii) one-particle 
one-hole (lplh) excitations, that is \ap) = e^e^HF), 
and (iii) two-particle two-hole (2p2h) excitations, i.e. 
\<jp, p\) = e'lej ) e fl e\\BF}, is used in order to diagonalize 
the Hamiltonian. Notice that a < p are single-particle 



states above the Fermi level, and p < A are states below 
the Fermi level. 

In the Hilbert subspace with the same quantum num- 
bers of the Hartree-Fock state, the electronic Hamilto- 
nian takes the form: 

/E HF D\ 
H = A B , (8) 
V -D* B l C I 

where E HF = (HF\H\HF) is the Hartree-Fock 
energy, j4 ct >',o> = ( (J '/ i 'l^l cr / i } is the Tamm- 
Dankoff matrix, DnF.appX = {HF\H\ap, /iA), 
Ba'^'^pnX = {<?' v'\H\ap,p\), and C a i p -^x\crpp\ = 
{a' p' , p'X'\B\ap, pX). D l and B l are, respectively, the 
transposes of matrices D and B. In sectors with quan- 
tum numbers others than the Hartree-Fock state, the 
first row and column of matrix © should be dropped. 
Explicit matrix elements are given in Appendix A of 
Ref. 0. 

The ground-state energy, E gs , in this case, is estimated 
as the lowest energy state in each of the computed 2p2h- 
intraband spectra. 

The dimension of the Hamiltonian matrix was con- 
trolled by using an energy cutoff in the excitation energy. 
The latter is expressed in terms of the confinement en- 
ergy. The energy spectra is obtained by exact diagonal- 
ization of matrices with dimensions of about 10 4 - 10 5 . 
Several computations of the energy spectra of different 
systems were carried out in order to check the conver- 
gence of E gs (and consequently, of E corr ) with the energy 
cutoff. 



C. Configuration Interaction Scheme (Exact 
Diagonalization) 

Using the eigenfunctions of the single particle problem, 
the Hamiltonian can be written in second quantized form 
as 

B ^ ^ ^ijC^Cjfj -\- —pint ^ ^ ^ijkm^ia-^jfj' ^rna' ^ka ^ 
i(J ijkmcra' 

(9) 

where en are the diagonal one-body energies, and Vijkm 
are the two-body Coulomb interaction elements. 

We consider the 120 lowest single particle eigenstatcs 
(15 oscillator shells), and truncate the many-body basis 
by allowing excitations only to configurations with lim- 
ited non-interacting energy, given by the first term in Eq. 
@ . The interacting Hamiltonian is then diagonalized by 
using the Lanczos algorithm. 

The CI method results to a converging ground state 
energy as a function of the basis energy cutoff. It per- 
forms comparably to VMC in the six-particle case prov- 
ing the validity of the VMC results for the energies. A 
sample of these results can be seen in Fig. Q]for the case 
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FIG. 1: (Color online) Covergence of CI energies as a function 
of the energy cutoff (lines are guides to the eye). The HF and 
VMC results for E ga are also shown. All results corresponds 
to a system with TV = 6 and Tiuj — 2.78 meV. 



of hui = 2.78 meV. When the number of particles is 
increased, the CI basis size grows exponentially. Thus, 
larger systems are studied by means of the VMC method. 



This leaves us with the two variational /3o- ii0 -, parame- 
ters, one for parallel and the other for antiparallel spins. 
We optimize these parameters by minimizing the vari- 
ational energy with the stochastic gradient approxima- 
tion, [m 



III. RESULTS AND DISCUSSION 

We show in Table Q] the HF and VMC energies of 
dots with N = 6,12,20,30,42 and 56, and confinement 
strengths that range from 2.78 to 1110 meV. All the sys- 
tems considered in this work are closed-shell QD's with 
ground-state angular momentum and spin quantum num- 
bers L = S = 0. 

Notice that, for the parameters used in the calcula- 
tions, the scaled variable z = V 1 / 4 /? takes values in the 
range < z < 5, i.e. the whole weak coupling and part 
of the strong coupling intervals, z w 1 is the transition 
point from the strong to the weak coupling regimes [24j . 
GaAs parameters, m = 0.067 mn and e = 12.8, were used 
for numerical computations (25[. 

In the following, VMC results, which have proven to 



D. Density Functional Theory (Local Density 
Approximation) 

Our density-functional calculations employ the local- 
density approximation parametrization by Attaccalite et 
al [I?} and a Bessel function basis [l8|. The system di- 
mension, grid spacing, and the number of basis functions 
are chosen carefully to guarantee numerical convergence 
of the energies. 



E. Variatonal Monte Carlo 

We use, in the variational Monte Carlo pjl] calcu- 
lations, wave functions of the Slater- Jastrow ty pe [20l l. 
that arc shown to be accurate for quantum dots. 



22] The 



Slater part, corresponding to a nonintcracting system, is 
the product of two determinants: one for spin-up, an a 
second for spin-down particles. The Jastrow factor is of 
the form 



exp[J] = exp 
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(10) 



are fixed 



where 

by the cusp conditions. In two dimensions, they are, 
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TABLE I: Hartree-Fock (E H f) and VMC (E VM c) energies 
(in meV), of the systems considered in this work. 



N 


Hio [meV] 


Ehf [meV] 


Evmc [meV] 




2.78 


82.0805 


77.6597(9) 




5.55 


136.319 


131.032(1) 


6 


11.1 


230.152 


224.067(2) 




111.0 


1519.82 


1511.82(2) 




1110 


12439.3 


12429.5(0) 




2.78 


270.841 


262.622(5) 




5.55 


446.772 


435.330(5) 


12 


11.1 


743.298 


730.172(1) 




111.0 


4622.29 


4604.75(9) 




1110 


36070.0 


36047.2(3) 




2.78 


653.023 


635.146(3) 




5.55 


1063.55 


1043.43(1) 


20 


11.1 


1755.28 


1732.28(2) 




111.0 


10522.4 


10491.3(2) 




1110 


79410.8 


79367.5(7) 




5.55 


2110.15 


2078.81(1) 




11.1 


3464.53 


3428.73(2) 


30 


22.2 


5768.21 


5728.49(3) 




33.3 


7833.70 


7791.35(2) 




111.0 


20244.4 


20195.1(3) 




5.55 


3719.60 


3674.32(4) 




11.1 


6084.43 


6032.74(1) 


42 


22.2 


10082.0 


10024.5(2) 




33.3 


13647.5 


13587.4(2) 




111.0 


34875.4 


34805.1(9) 




8.33 


8022.79 


7955.46(6) 




11.1 


9841.43 


9770.79(7) 


56 


22.2 


16247.6 


16169.2(7) 




33.3 


21937.6 


21853.9(7) 




111.0 


55558.3 


55459.5(1) 



N=6 N=12 N=20 N=30 N=42 N=56 
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FIG. 2: (Color online) Accuracy, E l , of the ground-state en- 
ergies calculation as a function of ft for the different compu- 
tational methods i used in this work. Lines are guides to the 
eye. 
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FIG. 3: (Color online) Correlation energies (in units of the 
confinement energy Tioj) as a function of ft for each particle 
number iV and used computational methods. Lines are guides 
to the eye. 



be very accurate [23j, were taken as references. We define 
the relative deviation, 
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The index i in Eq. (|12[) runs over the used methods. 
Fig. [5] shows E 1 as a function of j3 for each particle num- 
ber and for each method. 

From Fig. [2] it can be seen that, FCI practically co- 
incides with VMC for N = 6, whereas DFT and VMC 
results are very similar for larger N. The maximum con- 
tribution of E corr to E gs is less than 6%, and corresponds 
to the case of N — 6 particles with the weakest confine- 
ment. 

In Fig. |3l we show E corr in units of the confinement 
energy Hoj. The amount of correlation captured by the 
2p2h CI approach is significantly lower than the amount 
captured by other methods, in particular for systems with 
large TV under weak confinements. Let us stress the (ex- 
pected) increase of E corr with /3. 

In the following, we use VMC results in order to es- 
timate E corr because (i) the FCI scheme can not be ex- 
tended to systems with large particle numbers, and (ii) 
the DFT(LDA) is not a fully variational method. A uni- 
versal relation of the form E corr /Tiuj = N a / CO rr(TV 1 / 4 /3) 
is tried. We observe that results for different N and (3 
can be accomodated along a single curve when a « 3/4 
(Fig. 0]). For the function f corr , a two-parameter fit of 
the form az 1 leads to: 



The scaled correlations energies are shown in Fig. |4] 
along with the function / corr (z). It can be seen that the 
quality of the obtained analytical expression for f corr is 
very good in the whole range of z considered, showing a 
mean deviation of about 5% for E corr . 

Eq. p3[) can be straightforwardly rewritten in such a 
way that the direct dependence of E corr on the system 
parameters (say, N and Tilo) become explicit. For the 
specific case of GaAs it reads 

E corr (N,hu}) = -0.446 (huj) 1/e N 7/e , (14) 

where E corr is given in mcV. The correlation energy per 
particle shows exponents equal to 1/6 for both Tiui and 
N. 

Another interesting quantity is the fraction of the total 
energy corresponding to E corr , 



X = 



E, 



(15) 



It can be shown that x a l so follows a similar scaling 
relation as a function of the parameter z. We use Eq. 
(1131) in the numerator and the estimation 



E 



_ 2 0.698 z + 1.5 z 4/3 + 2.175 z 5 / 3 
huN 3 / 2 ~ 3 + 1 + 2.149 z 1 / 3 + 1.5 z 2 / 3 + 2.175 z ' 

n • (16) 

obtained in Ref. [3[, for the denominator. This process 
leads to 



= -0.060 z^ 3 . (13) X (z)N^ = f x (z) = ^ (17) 
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FIG. 4: (Color online) Scaled correlation energies as a func- 
tion of the variable z = (3N 1 ' 4 . The solid curve represents 
the function f co rr in Eq. (|13p . Some results from Ref. [2|| are 
also included. 



FIG. 5: (Color online) Scaled relative correlation energies as 
a function of the variable z = /3N 1 / 4 . Symbols are the results 
from the VMC calculations, and the solid line corresponds to 
Eq. (|17[) . Some results from Ref. [2(J are also included for 
comparison (crosses). The Inset amplifies the small z region. 



where 

p(z) = 0.18 z 5/3 + 0.387 z 2 + 0.27 z 7/3 + 0.392 z 8/3 (18) 
and 

q(z) = 2 + 4.298 z 1/3 + 3 z 2/3 + 6.444 z 

+ 4.5 z 4 / 3 + 6.525 z 5 / 3 . (19) 

The scaled computed values of \ are shown in Fig. [SJ 
along with the obtained analytic expression. It can be 
seen that Eq. (fl7|) works remarkably well for systems 
with large particle number. This is because the expres- 
sion used for the ground state energy, i.e., Eq. (|16|) . was 
originally obtained in the large- N limit. 

Notice that, in the strong confinement region (z < 1), 
Eq. (fl~7j) is able to describe all of the systems, no matter 
the particle number (See inset in Figj5]). On the other 
hand, in the weak confinement limit (z 3> 1), Eq. (fl7|) 
predicts a linear dependence of the scaled x on the pa- 
rameter z. From this linear dependence we can obtain a 
rough estimation of the behavior of \ (for z > 1) as a 
function of TV and Tuo, that is, x( z ^ 1) oc 1/ \/N(huj). 

The quantity \ can be used to study the role of corre- 
lations in artificial atoms. It is known that they become 
relevant for systems with small particle number and weak 
confinements. Our calculations support this known re- 
sult. However, even for any large-iV system, Eq. (jTTJ) 
shows that we can find a small enough value of hiu lead- 
ing to a given degree of correlation. 

We finally consider some previous calculations of both, 
E corr and \ m order to establish comparisons. For in- 
stance, in Ref. [26[ computed values of these quan- 



tities for QD's with TV = 2,6, 12 and 20, and con- 
finement potentials between 3.11 and 11.1 meV were 
reported. These results were obtained from different 
Coupled-Cluster methods. Our results are in agreement 
with those in Ref. [26|, as it can be seen in Fig. [4] and 
Fig. El 



IV. CONCLUDING REMARKS 

In summary, we have performed extensive numerical 
calculations for charged semiconductor quantum dots 
and found a universal scaling relation for the correlation 
energy, Eq. ([3]) , which resembles the scaling coming from 
TF theories. This is an unexpected and very interesting 
result which, in principle, could have direct implications 
in the design of new correlation functionals for DFT cal- 
culations [27j . 

A universal scaling relation for the fraction of the total 
energy associated to the correlations was also obtained, 
Eq. (|17p . Such expression would provide information on 
the degree of correlation of the system and the accuracy 
of the HF estimation, even without any calculation. The 
material parameters (effective mass, dielectric constant) 
are contained in the scaling variable z. The numerical 
coefficients entering the r.h.s. of this equation, on the 
other hand, do not depend on the material the system is 
made of. 

We notice that in real atoms, where TF theory predicts 
for the total energy the dependence [28| : 

E gs (N,Z)^N J / 3 f gs (N/Z), (20) 
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where Z is the nuclear charge, the correlation energy 
seems also to show a scaling a la TF: 

E corr (N, Z) « N a f corr (N/Z). (21) 

The coefficient a is near 4/3 [29U38j ] . We stress that real 
atoms, unlike artificial ones, are always under a weak 
correlation regime. Research along this directions is in 
progress. 
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